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Abstract 

We show that a form for the second partial derivative of 1/r proposed by Frahmi and subse- 
quently used by other workers applies only when averaged over smooth functions. We use dyadic 
notation to derive a more general form without that restriction. 
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In Ref . HI the equation 

9id ' i;J = " [f ) M(r) + '=^pr^ (i) 

is proposed for the second partial derivative of 1/r. This result has subsequently been used 
in more recent papers^ - -. The purpose of this comment is to show that the derivation of 
Eq. ([I]) in Ref. 1 is flawed, and to present a direct derivation of this second partial derivative. 
While Ref. 1 uses two indirect methods to deduce Eq. flTJ, we employ dyadic notation^ to 
take the partial derivatives directly. 

In dyadic notation, the left-hand side of Eq. ([I]) can be operated on as 

W (I) = -V (L) (2a) 
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where I is the unit dyadic. 

To evaluate the term V(l/r 3 ), we start with the well known identity 



V • 



(£)=4**(r). (3) 



The left-hand side of Eq. (j3J) contains no vector other than r. Therefore the scalar functions 
on each side of the equation can have no angular dependence. The identity of Eq. (jHJ) is 
usually proven by integrating the left-hand side over a volume, and applying the divergence 
theorem. The volume integral over the left-hand side becomes an integral over the bounding 
surface. Since the integrand has no angular dependence, the integral over the solid angle 
equals 4-7T if the volume contains the origin. It equals zero if the volume does not contain 
the origin. Thus, the right-hand side of Eq. (jHJ) equals 4n5(r) by the definition of the delta 
function, and the identity is proven. 

We can express the left-hand side of Eq. (j3J) as 

V-f^^+r.vfiUi + r.vfiy (4, 
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which isolates the term V(l/r 3 ). Since the function (1/r 3 ) depends only on r, its gradient 
must be in the f direction. Thus we can write 

V Qj) = ^(r), (5) 



where the scalar function g(r) is give by 



y(r) = v V(l). (G) 



We combine Eq. (JBJ) with Eqs. ([3]) and (BJ to find 



and then 
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Substituting Eq. ((8]) into Eq. ( |2cl) . we obtain 



VV ( - ) = ^-\-Aizrr5(r). (9) 



In the Cartesian tensor notation of Ref. 1, Eq. ([9]) would be written as 

\ f J f 

which differs from Eq. ([T]) in the delta function term. 

Objections have been raised about the relevance of Xi/r multiplying a delta function,- 
because Xj/r is not well defined in the limit r — > 0. However the same objection could be 
raised against the delta function itself, which is also undefined at the origin. As with the 
delta function, factors of x^/r give a definite result when used in a volume integral, even 
when multiplied by a delta function. Also, our Eq. ((8]) shows that the gradient of 1/r 3 could 
not be written without r (with Cartesian component Xj/r) multiplying the delta function. 
The ratio Xi/r is used in several places in Ref. 1 for infinitesimal r. It is inconsistent to 
preclude x^/r in one place and then use it in another. 

Where does Ref. 1 lose the Xj/r factors multiplying the delta function? Reference 1 
arrives at its expression for didj(l/r) by first using a "plausibility argument" and then a 
"physicist's proof." The plausibility argument depends on the statement: "Noting that the 
mixed second derivatives cannot contain a delta function. . . ." The falsity of this statement 
is demonstrated by our result in Eq. (TTUT) which shows that the mixed second derivative does 
contain a delta function, and does not depend on how the first derivative was taken. The 
plausibility argument arrives at the form 5« instead of XiXj/r 2 by incorrectly assuming it. 
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The physicist's proof in Ref. 1 uses an integral over solid angle to deduce that the form 
in Eq. ([T|) is correct. This proof uses the identity 

for the integral over the surface of a sphere of infinitesimal radius R. In dyadic notation, 
this identity is 

ffdQ = — 1 (12) 
3 v ; 

Reference 1 then uses Eq. ( ITT)) to show that its form for the delta function term gives a 
correct integral when multiplied by "an arbitrary smooth function" and integrated over a 
sphere with infinitesimal radius. This derivation works because a function that is smooth at 
the origin has a Taylor expansion which has no angular dependence in the limit of vanishing 
radius, and because Eq. (ITT)) shows that our form for the delta function term reduces to 
that of Ref. 1 when averaged over solid angle. However this averaging procedure restricts 
the applicability of the form in Ref. 1 to smooth functions. For instance, a function such as 
f(p, f) = (p- r)(p ■ f), which is not smooth at the origin, gives a different result in a volume 
integral when multiplied by 15 (r) instead of rr5(r). That is, 

f 47T S 87T 

(p-r)(p-r)rr5(r)dV = —1+ —ppi (13) 

while 

f a 47TJ 

(p-f)(p-r)1 5{r)dV = —I. (14) 

Even though it is incorrect, the use of Eq. (1) in most physics applications leads to the 
correct result, since functions used in physics are usually smooth at the origin. For instance, 
most electromagnetism textbooks derive the singular part of the electric field of an electric 
dipole p by applying the divergence theorem, and effectively averaging E over all solid angle. 
In this case, either equation ([1]) or (Q gives the same result (— A7ip/3)5(r) for the singular 
part of E. However, writing 



E = 3 (r-')'-P _ *ZR 5{r) (15) 



is inconsistent mathematically, because the delta function term is averaged over solid angle, 
while the first term is not. The use of Eq. ([9]) leads to the mathematically consistent 
equation^ 

E = ^PjJ^lJP -4n(p-r)r5(r). (16) 



The magnetic field of a magnetic dipole /x also involves the gradient of 1/r 3 as can be 
seen by writing 8 

B = VxA = Vx O^r) (17a) 

= ^ ■ (£) " » ■ V (£) ( 17b ) 
= MV.(^)-^-r M -v(i) (17c) 
3(/Lt • f)r — /x 



47r(/x • f)r5(r) + 47rjx5(r). (17d) 



For the last step we have used Eq. for the divergence of r/r 3 , and Eq. (JSJ) for the gradient 
of 1/r 3 . We see that the field of a magnetic dipole is like that of an electric dipole, but has 
an additional singular term 47r/x5(r). The singular part of Eq. ( 117dl) . averaged over all solid 
angle, is given by (+87r/x/3)5(r), which is the form given in most textbooks for the singular 
part of the magnetic field of a magnetic dipole. 

In summary, we have shown that the second partial derivative of 1/r can be found by 
direct differentiation using dyadic notation. As a mathematical statement, we have shown 
that Eq. (1), as proposed by Ref. 1, cannot be used with functions that are not smooth at 
the origin. The higher derivatives considered in Ref. 1 and Refs. |2|45| would also be affected 
by using Eq. rather than Eq. (1). (We do not consider higher derivatives here to keep 
the paper relatively simple and accessible.) 
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